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ABSTRACT 

A model for Type II solar radio outbursts is discussed. The radiation is 
hypothesized to originate in the process of collective bremsstrahlung emitted by 
non -thermal electron and ion density fluctuations in a plasma containing a flux of 
energetic electrons. The source of radiation is assumed to be a collisionless 
plasma shock wave rising through the solar corona. It is also suggested that the 
collisionless bow shocks of the earth and planets in the solar wind might be 
similar but weaker sources of low frequency emission with a similar two- 
harmonic structure to Type II emission. 
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AND TYPE H SOLAR OUTBURSTS 


by 
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INTRODUCTION 

In a previous paper 1 we calculated the bremsstrahlung emitted from thermal 
plasmas which co-exist with a flux of energetic (suprathermal) electrons. It was 
found that under some circumstances the radiation emitted from such plasmas 
can be greatly increased compared to the emission from a Maxwellian plasma 
with no energetic particles present. The enhanced emission occurs at the funda- 
mental and second harmonic of the electron plasma frequency. 

The object of this paper is to discuss a model for Type II solar radio out- 
bursts in terms of the above radiation process. It has been suggested by several 
authors (for a review of the literature see reference 2) that this radiation is 
probably produced by longitudinal plasma oscillations which are excited by some 
mechanism associated with a shock wave rising through the solar corona. The 
longitudinal waves then convert part of their energy 3,4,5 into electromagnetic 
radiation with a frequency co e by scattering on density fluctuations in the corona, 
and part into radiation at 2 co e by scattering of longitudinal waves by each other 
( a) e is the electron plasma frequency (47m 0 e 2 /m) 1/2 ). 

The present model follows essentially along these same lines with the new 
feature that we are able to calculate quantitatively the bremsstrahlung from an 
assumed distribution function for the suprathermal electrons in the shocked 
plasma and point out the importance of Landau damping of the longitudinal waves 
by the suprathermal electrons. In the absence of adequate theories of collision- 
less plasma shock waves to guide us, we also assume that the shock wave in the 
corona has some of the structural features experimentally observed to charac- 
terize the earth's bow shock wave in the solar wind. 6 

One can visualize the energetic electrons behind the shock front as exciting 
Cerenkov electron plasma waves which then, through scattering from each other 
and by ion waves, convert part of their energy into electromagnetic radiation. 
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We demonstrated previously, 1 however, that the emission intensity is very sen- 
sitive to the velocity dependence of the distribution function, f (v) , in the high 
energy tail (f includes both thermal and suprathermal electrons). This is because 
the Cerenkov electron plasma oscillations are also re-absorbed through Landau 
damping by the energetic electrons as well as emitted by them. This leads us to 
suggest a distribution function of a certain type for the energetic electrons in the 
shock wave (see Figure 3), namely one in which there is a gap in velocities be- 
tween the thermal and energetic electrons. Such distributions are highly radiative 
(Note since f = f( |v |) they are also stable) . In contrast, distributions which 
have high energy tails on them (representing a flux of energetic electrons in the 
medium) which merge continuously with the "thermal" particles, emit very little 
electromagnetic radiation. An example of this, namely f e = 4V 3 /(v 2 + V 2 ) 3 77 2 , 
was treated in reference 1. 

We shall next summarize briefly some of the observed properties of Type II 
outbursts. It is then also suggested that the earth's collisionless bow shock might 
emit electromagnetic radiation with a similar structure to Type II solar emission, 
but at lower frequencies (f = co /2 tt - 20-300 kc./sec.). This radiation would 

e 

only be detectable outside the magnetopause. 


CHARACTERISTICS OF TYPE II OUTBURSTS 

The data in this and the following section is taken from several sources. 

They include a review of solar radio emission by Roberts, 7 the chapter by Pawsey 
and Smerd 8 in "The Sun," papers by Maxwell and Thompson, 9 and more recent 
reviews of the literature by Wild, Smerd, and Weiss, 2 and by Boyd. 10 

Type n, or "slow -drift" radio bursts, appear on the average about 6 minutes 
after the commencement of a solar flare and probably originate above about 
10 s km. over the photosphere. They are unpolarized and a large outburst 9 radi- 
ates a total energy integrated over frequency and time of order 10 22 ergs. This 
is as if an area of 10 20 cm 2 radiated like a black body at a temperature of 10 13 , 
i.e., they are a factor of about 10 7 times more intense than the normal solar 
thermal emission from the above area. 

The duration of these bursts is approximately 10 minutes, and their fre- 
quencies drift from high to low frequencies, which is interpreted as due to the 
radio source moving up through the solar corona at about 1 to 1.5 10 8 cm/sec. 
Their velocity is deduced from the frequency drift by assuming that the fre- 
quencies observed correspond to the local plasma frequency at the source, which 
in turn depends on knowing the electron density in the region. 
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Independent measurements of the height of the radio source in the corona 
using a swept-frequency interferometer 2 leads to electron densities, n Q , at the 
source which are of order 10 times the accepted values for the quiet corona. 

The deduced variation of n Q with height agrees closely with those deduced opti- 
cally by Newkirk 11 for typical coronal streamers. The source is thus probably 
propagating out along a coronal streamer. (See Figure 1.) 

One of the most puzzling features of these radio outbursts is the existence 
in most cases (80%) of a strong second harmonic in the spectrum. Occasionally 
the second harmonic may be even stronger than the first, and yet there is no 
evidence of higher harmonics, even though they would have been visible in Maxwell 
and Thompson's work 9 at an intensity of 20 dB lower than the fundamental or 
second harmonic. The maximum starting frequency, f , of the fundamental rarely 
exceeds 175 Me. and often drifts to below 25 Me. (It should be noted that most 
experimental authors quote frequencies in terms of the - angu l a r frequency f = co e /2rr) . 

The two harmonics sometimes have a similar structure, i.e., variations in 
intensity or line breadth occur simultaneously in both harmonics. They also 
sometimes both split into two sub-bands with a splitting frequency small com- 
pared to co e . This latter effect is generally thought to have its origin in coronal 



Height (km.) above photosphere. 


Figure 1-Plot taken from Maxwell and Thompson (Ref. 9) of 
plasma frequency and electron density as functions of height 
above the photosphere. The lower curve is based on the 
Baumbach-Allen (B.A.) density distribution for the undisturbed 
corona, and the upper curve (1 0 x B.A.) are values for a corona I 
streamer along which a Type II burst propagates. 
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magnetic fields. Occasionally a "herring bone" effect is also observed, 12 i.e., 
small fast drift bursts (Type Ill’s) accompany each of the harmonics. 

Another feature recently observed in some events is that the second harmonic 
appears to originate from deeper in the corona than the fundamental. 13 Smerd, 
Wild, and Sheridan 14 suggest that this is due to refraction of that part of the 
second harmonic radiation which is emitted down into deeper regions of the solar 
atmosphere. Additional refraction effects which may also play a role are due to 
the radiation traversing the irregular shock front and in leaving the coronal 
streamer and entering the quiescent corona. 


MODEL FOR TYPE H OUTBURSTS 

Suppose we assume that in the course of a solar flare a shock wave is gen- 
erated by the upward expansion of the hot flare plasma into the solar corona. 

In the absence of adequate theories to guide us we shall assume that the collision- 
less shock wave is similar in structure to the bow shock wave which has been 
experimentally observed to be generated by the obstruction the earth's magnetic 
field presents to the "super-Alfvenic" solar wind. 6 Thus the shock transition 
is likely to be turbulent in nature and to propagate at some arbitrary angle to 
any undisturbed relatively ordered coronal magnetic field which may exist ahead 
of the wave (see Figure 2). We note that the shock velocity (of order 10® cm/sec) 
is considerably greater than the Alfven speed, B // 477Mn Q , for the range of den- 
sities and field strengths involved at the source of emission in the corona. 

(B = 0(1 gauss) and n g 'k 10 6 - 10 8 cm -3 )* 

Now in the turbulent region behind the collisionless bow shock wave of the 
earth, fluxes of energetic electrons are observed. 15 The mechanism by which 
they are energized in or behind the shock front is not yet clear, although it is 
most plausible that it involves a stochastic process 16 in the plasma turbulence 
behind the front of the wave. We assume that in the collisionless shock propag- 
ating up through the coronal plasma, a similar condition of disordered plasma 
and magnetic field co-existing with a flux of suprathermal electrons exists. On 
the basis of these assumptions we shall apply our bremsstrahlung calculations 
to see how many of the spectral characteristics of Type II's we can explain. 


*The magnetic field strength in the outer corona is not well known. In reference (9) it is pointed 
out that if the band splitting in Type II emission is interpreted as a gyromagnetic effect analogous 
to Zeeman splitting, the following field strengths are obtained: 1.5 gauss at 8.10 5 km. above the 
photosphere, 3 gauss at 4.10 5 km., and 8.5 gauss at 3.10 5 km. 
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Figure 2-ldealized sketch ofa model for a collisionless shock wave inthe coronal plasma 
above a solar flare (assuming similar characteristics to the earth’s bow shock) 


The turbulent magnetic field behind the shock front will serve both to con- 
fine the energetic electrons to that region, and through scattering constrain the 
velocity distribution for the suprathermal electrons to be approximately isotropic. 
Although some energetic electrons are likely to escape into the ordered field 
ahead of the shock, they will rapidly travel upwards along the field lines and 
thus their intensity ahead of the shock should be much less than that behind the 
shock front. 
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Previously we calculated the collective bremsstrahlung 1 emitted from 
plasmas containing energetic particle fluxes, using several distribution functions. 
For the present application we shall assume that the plasma behind the shock 
front in the rest frame of the shock can be represented by the following distribu- 
tions, f. and f for the ions and electrons, 

l e 


f . (v) = exp 

(277 v ?) 3/2 



( 1 ) 


f (y) = £ exp (--Zl-V (1 -P) I s (ly| - V E ) g(|y| ), (2) 

(2 ,V e f /2 ^ 2V g 2 J 

where I = 0 if I v I < V_ and = 1 otherwise. V and V. are the electron and 
ion thermal velocities, V = (KT/m) 1/2 = (M/m) 1/2 V. . We also choose 
exp(-V 2 /2V 2 ) « (1 - /S) « yS ~ 1 and V E > several- V e (see Figure 3). Note 
the distributions are all normalized to unity, J~J^ dv f e = 1, and the number 
density of the electrons and ions is n 0 . The function g( |v | ) represents the 
distribution for the suprathermal electrons. 

Now f e (v) as given in (2) represents a Maxwellian gas of electrons co- 
existing with a tenuous isotropic flux of energetic electrons all of which have 
velocities |v | > V E . We found that distributions of this type are highly radiative 1 
and give enhanced emission at a> e = (47m 0 e 2 /m) 1/2 and 2a> e . The basic reason 
for this is due to the gap in velocities V < | v | < V E which only the few Max- 
wellian tail particles occupy. Electrostatic (longitudinal) Cerenkov waves are 
emitted into the phase velocity range several- V g < co/k < V E by the energetic 
electrons, but only reabsorbed (Landau damped) by the few thermal particles 
far out on the Maxwellian tail. These longitudinal waves then scatter off of each 
other and emit electromagnetic radiation. 

A distribution function of this type is expected if the acceleration mechanism 
that energizes a preferred few of the electrons in the plasma is very efficient in 
a velocity range | v | < V , so that they rapidly get carried up into the range 
|v|>V E . 
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Figure 3-Schematic plot of f ( | y |)as a function of |v_ |for/S < 1 as given by Equation (2) 

On the other hand a resonance distribution, f . ( I v I ) = 4V 3 A/fv 2 + V 2 ) 3 , 

for which the suprathermal electrons do contribute appreciably to Landau damp- 
ing for all longitudinal waves of any phase velocity, gives little increase in the 
emission. 1 This feature is characteristic of distributions for which the distri- 
bution of the "suprathermal" particles merges continuously in velocities with 
that of the "thermal" particles. 

For the distributions (1) and (2) the first and second harmonic emission is 
(taking C = 1/2 in eq. (58) ref. 1). 


n ff\ /mV 2 \ /V \ 3 i / V \ 

1 (ipo - epo) = zEi G J [ — ] <w e(“) — ( 

277/3 \n 0 L 3 y e \ c j V\\j 


(3) 


and 


e rgs/sec./cm 3 . 


I a „ (epo-epo )^ 771 ^^ 1 zJUL G 2 J 2 


m V 2 \ /V \3 . /V \ 2 

e \ j e \ 1 

n Q L 3 y 6 \ c / L 3 \ 'Ey 


(4) 
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where, 


Ji 



exp 



(5) 


J 2 




( 6 ) 


and 


G = 



V E v dv g(v) , 


( 7 ) 


and L is the Debye radius (L = V e /a> e ) . Note that the magnetic field has been 
neglected in calculating the above results. However although the coronal mag- 
netic fields play a dominant role in the shock structure, and probably also in 
energizing the electrons represented by g( |v |), they are of secondary importance 
in computing the bremsstrahlung at frequencies a > = 0(co e ) >> Q = eB/mc . We 
are making calculations to see if the band splitting can be accounted for by in- 
clusion of a "weak" magnetic field. 

The emission intensities (3) and (4) were derived from convolutions involving 
the spectral densities S (k , aj) and S. . (k , co) for the fluctuating number densi- 
ties of the electrons and ions. For wave numbers |k| < k D = L -1 (k D is the 
Debye wave number), the spectral density S ee has resonances at co = 0 , a> e , 
and s. . has a resonance at co = 0. These resonances (see Figure 4) correspond 
to a spectrum of electron plasma oscillations (epo's) and ion plasma oscillations 
(ipo's) in the medium in the sense that they represent Fourier components of the 
spectral densities for which co and k are approximately related as they would be 
for propagating Vlasov ion and electron plasma waves. Thus, roughly speaking, 

I derives from epo-ipo scattering, and I from epo-epo scattering. 

u^ e 2 co 
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Now we note (Figure 4) that the presence of energetic electrons in the plasma 
enhances the epo components of S ee at ^ = o> e , but the fast electrons A do not pro- 
duce a corresponding excitation of the low frequency ipo modes at o> = 0 . (Com- 
pare S for y8 = 1, J3 < 1 in Figure 4). This results in a situation in which if 
the exponential factors exp(V^/2V^) in (5) and (6) are large, then l 2a) » I a . 

This is contrary to the observed fact that 1^ and 1 2 <a e are of comparable orders 
of magnitude for a Type II outburst. In order to explain the observations we must 
find a mechanism which does enhance the low frequency mode excitations in the 
plasma, and the corresponding I . 



Figure 4-Schematic plot of the spectral densities S ee and S n as functions of frequency 

for a fixed k < kp 
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The most obvious source of low frequency density fluctuation is the shock 
turbulence itself. Now we are not yet in a position to calculate the spectral den- 
sity for this turbulence - what is needed is a theory of turbulence for Vlasov 
plasmas and shock wave structure. However in the case of the earth's bow shock 
the magnetic field fluctuations are observed to be 6 such that | AB | is of 0( | B| ), 
and similarly | 8/o| = 0(1) for the normalized density fluctuations and for wave- 
lengths greater than an ion cyclotron radius. In the theory of reference 1, how- 
ever, the thermal fluctuations of a "homogeneous" plasma are such that | hp\ = 0(g) 
where g = l/n Q L 3 is the plasma parameter. For the coronal plasma g ~ 10 -8 . 
Thus if the coronal shock wave is similar in structure to the earth's bow shock, 
then the spectral densities at frequencies & << co will be increased by many 
orders of magnitude, perhaps by as much as 10 8 in some wave numbers over 
those used in the derivation of (3). This suggests that it is capable of bringing 
I up to the intensity level of I , but of course does not enable us to make a 
quantitative calculation. Since in our model I 2o) most directly reflects the 
simpler properties of energetic electron fluxes in the shocked plasma, we shall 
in the rest of this section, examine the conditions under which I 2a>e is large 
enough to account for the observed power level at 2^ e . 

Consider for example the parameters (see Figure 1), T = 10 6 °K, V e /c = 10 -2 , 
n Q = 10 8 , f ~ 100 Mc/sec, i.e., w e = 600 Mc/sec., and for the Debye length, 

L = V /cj - .5 cm. Using these in (4) we find, 

e e 

I 2ct) = 2.10 -14 (1 - /3) 2 G 2 exp e rgs/sec. /cm 3 . (8) 


Various assumptions are now possible for the flux of suprathermal electrons 
behind the shock front. Thus, suppose we assume that the energy density of the 
suprathermal electrons is comparable to the thermal energy density of the plasma, 
i.e., 


(1-/3) n Q m 


= Z 3 n 0 m v e 


Further, from the observed typical linear dimensions of solar flare plasma, we 
take the radiating volume of plasma behind the collisionless shock front to be 
~ 10 29 cm 3 . Equation (8) now gives for the total emission (note G ~ 0(1/477), /3~1), 


jtotal ~ -i q13 

2^ e 




ergs/sec. 


(9) 
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If the energetic electrons all lie in a velocity range |v| > V E = 5.5 V e (i.e., are 
in the energy range > 2.5 keV), the total emission becomes I £° tal = 2. 10 19 
ergs/sec., which is the required power output to explain type II outbursts. 

In reference 1 we noted two conditions for the applicability of (4) and we 
shall now show that these are satisfied. 


Landau Damping 


The first condition is that there be an equilibrium established between the 
emission and re-absorption (Landau damping) of longitudinal (electrostatic) waves 
for the small wave-number part of the spectral density in the radiating volume 
of plasma. This is the case if the linear dimensions, L s , of the shocked plasma 
is much larger than the propagation length (co e / k) y- 1 (k) for a wave of phase 
velocity a> e /k = V £ . (y L is the Landau damping decrement) . We require, 

= 3 km. (10) 



for the parameters we have chosen. This condition is therefore satisfied. 


Collision Damping 


The second condition is that Landau damping be more important than colli- 
sion (correlation) damping for electron plasma oscillations in the wave-number 
range considered. Equation (61b) of reference 1 gives this condition as, 


> J_ in (n 0 L 3 ) (11) 

n o L3 

For the parameters we have chosen the left side of this equation is about 10 
times the right side so that this condition is also met. We should point out in 
passing, however, that even if it were not satisfied, one would still expect greatly 
enhanced emission although our simple formula (4) would become modified. 



11 


LOW FREQUENCY EMISSION FROM THE EARTH'S BOW SHOCK 


For the purpose of estimating the emission in this case 6 we assume the 
following parameters: L 3 = 10 29 cm 3 = radiating volume, n Q = 10, T = 10 s , 
and (1 - (3 ) n Q m V 3 = a/3 mV 2 n Q , where a is approximately the fraction of the 
thermal energy represented by the energy density of the suprathermal electrons. 
Thus V = 1.2 10 8 cm/sec, co e ~ 1.5 10 s andL ~ 10 3 cm. Again we assume 
the electron acceleration mechanism carries a few thermal electrons rapidly up 
to velocities | v | - V £ , and that the region | v | < V £ is occupied only by a few 
Maxwellian electrons. In the previous section we made a similar assumption 
for the plasma behind a coronal shock wave since only then could we get sufficient 
emission to explain the observed intensity of a Type H burst. In the case of the 
region behind the earth's bow shock a similar energy gap between the thermal 
and suprathermal electrons may or may not exist, but we shall explore the hypo- 
thesis that it does. 

The second harmonic emission then follows from (4) as, 


I‘° tal Y 2.10* 6 a 2 

e 




ergs/sec. 


This is now extremely sensitive to the ratio (V e /V E ). Thus if for example we take 
V E /V e ~ 5.5 as before, and assume a = 0(1), we find I^ 31 = 500 ergs/sec. (at 
a steady frequency of f =£ w e / 2tt = $0 kc./sec.). The two e conditions (10) and (11) 
for the applicability of equation (4) are also satisfied for this choice. 

Such radiation if it exists would have to be observed outside the magnetopause. 
Although its power level might be quite low and one would have to contend with 
possible low frequency background noise from the Sun and magnetosphere of the 
earth, we think it is worth looking for such emission - both from the earth and 
other planets. We cannot of course make a reliable prediction of the expected 
intensity for such a process due to the sensitivity of I 2a)e to the ratio V E /V e . 
However, such a pheonmenon if detected would give us a much closer look at the 
physical mechanisms involved in the corresponding solar Type II events. 
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